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We propose a novel mechanism enabling us to have a continuum bump as a signature of gamma- 
ray excess in indirect detection experiments of dark matter (DM), postulating a generic dark sector 
having (at least) two DM candidates. With the assumption of non-zero mass gap between the two 
DM candidates, the heavier one directly communicates to the partner of the lighter one. Such a 
partner then decays into a lighter DM particle along with an “axion-like” particle (ALP) or dark 
“pion”, which subsequently decays into a pair of photons, via a more-than-one step cascade decay 
process. Since the cascade is initiated by the dark partner obtaining a non-trivial fixed boost factor, 
a continuum 7 -ray energy spectrum naturally arises even with a particle directly decaying into 
two photons. We apply the main idea to the energy spectrum of the GeV 7 -rays from around the 
Galactic Genter (GC), and find that the relevant observational data is well-reproduced by the theory 
expectation predicted by the proposed mechanism. Remarkably, the relevant energy spectrum has 
a robust peak at half the mass of the ALP or dark pion, as opposed to popular DM models directly 
annihilating to Standard Model particles where physical interpretations of the energy peak are not 
manifest. Our data analysis reports substantially improved fits, compared to those annihilating DM 
models, and ~ 900 MeV mass of the ALP or dark pion. 


I. INTRODUCTION 

There is astrophysical and cosmological evidence that 
DM exists in the Universe (see, for example. Ref. [1]). 
Relevant observations, mostly rooted in its gravitational 
effects, can be explained by postulating new stable par¬ 
ticles, not belonging to particle species in the Standard 
Model (SM). With this situation, there is a tremendous 
amount of effort to detect DM candidates: 1) direct de¬ 
tection experiments by measuring recoil energy of nuclei 
scattered off by DM, 2) indirect detection experiments 
by observing signals stemming from DM annihilation or 
decay, and 3) collider searches by actively producing DM 
particles and observing associated collider signatures. 
Among those experimental efforts, satellite-based cosmic- 
ray detection experiments such as PAMELA [2, 3], AMS- 
02 [4, o], and Fermi-LAT [6, 7] have received particular 
attention due to their great sensitivity to cosmic-ray sig¬ 
nals, giving rise to better chance to have not only con¬ 
firmation of the existence of DM but the information for 
deducing DM properties. 

The Fermi-LAT collaboration has provided the public 
data based on their observations, and a 7-ray excess at 
O(GeV) coming from the GC has been found. In partic¬ 
ular, it has recently reported in Ref. [b] that the excess 
exists, even assuming different foreground/background 
models. The relevant program was initiated by Ref. [9], 
and their intriguing observation has been strengthened 
by a series of their follow-up analyses and other indepen¬ 
dent groups [10-20]. Unlike other photon excesses such 
as 3.5 keV line [21, 22], 511 keV line [23], and 130 GeV 
line [24, 25], this is characterized by a continuum bump. 
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The basic claim is that the 7-ray excess spectrum is suf¬ 
ficiently consistent with the expected emission spectrum 
from charged particles in the SM into which DM parti¬ 
cles are annihilated. More specifically, the GeV excess 
is well-accommodated by a DM scenario where a pair of 
DM particles with a mass of ^ 30 — 40 GeV annihilate 
into a bb pair with an annihilation cross section of {av) ~ 
2 X 10“^®cm^/s [17, 20]. As an alternative annihilation 
channel, lepton pairs have been studied as well in Ref. [18] 
where they pointed out the significance of the contribu¬ 
tions coming from the diffuse photons from primary and 
secondary electrons that are produced in DM annihila¬ 
tion processes. They further analyzed the data including 
the inverse Compton scattering and bremsstrahlung con¬ 
tributions from electrons, and found that the data is well- 
described by ~ 10 GeV DM annihilating into a U pair, for 
which the associated annihilation cross section is given 
by (av) « (1 — 2) X 10“^®cm^/s [18]. Moreover, Ref. [20] 
showed that the GeV excess can be reproduced by other 
heavy SM final states such as W'^W~iZZjhhItt with 
a DM particle of tudm « 80 — 200 GeV and (av) « 
(2—8) X 10“^®cm^/s, depending on the final state with the 
systematic uncertainties in the gamma-ray background 
modeling taken into account. In Ref. [20], it was also 
shown that gg/W^W~/ZZ/hh/tt final states can pro¬ 
vide a good fit to the excess with « 40 — 200 GeV 
and (av) « (1 — 8) x 10“^®cm^/s. 

We remark that some of the realistic DM models have 
been proposed and studied: for example. Refs. [27-32] for 
the bb final state through a Higgs portal type interaction. 
Refs. [33, 34] for final state. Refs, [ io-37] for DM anni¬ 
hilating to a pair of on-shell particles that subsequently 
decay into SM // pairs. Ref. [38] even for 2” pairs of 
SM // final states from on-shell mediator pairs through 
multi-step cascades, and Ref. [39] for generic model con¬ 
straints. Although the recent report from the AMS-02 
collaboration [40] has started to rule out the qq final 
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FIG. 1: The example “dark” cascade scenario under consid¬ 
eration. 


state dominant DM scenarios explaining the measured 
relic abundance [41], it is straightforward to invoke hy¬ 
brid scenarios where bb and £i modes are mixed together. 

One should notice that the astrophysical uncertainty 
in 7 -rays coming from the GC in conjunction with the 
background modeling for the emission in the inner galaxy 
is still large. In addition, pions from the collision be¬ 
tween cosmic-rays and gas [10-12, 14] and millisecond 
pulsars [10-12, 14, 16, 42] can be sources to the GeV 
scale 7 -rays. Therefore, they have been proposed as a 
different approach to interpret the excess although the 
relevant spectral shape appears too soft at the sub-GeV 
energy regime to accommodate the observed energy spec¬ 
trum [43]. When it comes to the morphological feature 
for the observed excess, it is extended to more than ~ 10 ° 
from the GC beyond the boundary of the central stel¬ 
lar cluster that could contain a large number of millisec¬ 
ond pulsars [18], and observed distributions of gas seem 
to give a poor fit to the spatial distribution of the sig¬ 
nal [18, 44, 4-5]. We finally point out that very recently, 
another non-DM interpretation has been suggested by 
Ref. [46, 47]. They basically came up with a new method 
to characterize unresolved point sources based on which 
the excess can be explained by a population of unresolved 
point sources, giving a distribution consistent with the 
observed GeV 7 -ray excess in the relevant region. 

Notwithstanding those potential issues, we here pro¬ 
pose another novel mechanism to attain continuum en¬ 
ergy spectra of the 7 -ray excess, positing the DM inter¬ 
pretation. We first remark that there could exist multi¬ 
ple DM species, and the DM models based on such a DM 
framework can give rise to not only non-trivial cosmologi¬ 
cal implications (e.g., ‘‘^assisted freeze-ouf [48]) but inter¬ 
esting phenomenology (see, e.g., ^^boosted DM' [49-51]). 
In this context, we assume a DM partner having a non¬ 
trivial and fixed boost, which could be achieved by the 
annihilation of another (heavier) DM. The DM partner 
is further assumed to undergo a > 2 step “dark” cascade 
decay, and emit, in the final step, a (lighter) DM particle 
along with an ALP or dark pion that subsequently decays 
into two photons. FIG. 1 schematically depicts the ex¬ 
ample “dark” cascade scenario of our interest. We shall 
discuss the minimality of this DM scenario in conjunction 
with the elaboration of formalism later. We emphasize 
that the relevant ALP or dark pion comes with a non¬ 


trivial boost distribution, typically rise-and-fall-shaped. 
As a consequence, the resulting photon energy spectrum 
becomes a broad distribution.^ 

A novel feature of this type of energy distribution is 
that the peak of the photon energy distribution is ex¬ 
actly the same as half the mass of the ALP or dark 
pion [52, 53]. In other words, such a peak position is 
robustly connected to a physical property. This is not 
the case for other interpretations such as the DM models 
directly annihilating to SM particles because their energy 
peak highly depends on models of parton showering, dif¬ 
fusion mechanism, and so on from which the final photon 
spectrum is generated. We emphasize that the proposed 
strategy is completely generic to be applicable to any 
continuum bump in cosmic-ray energy spectra, not re¬ 
stricted to the 7 -ray excess at hand even if we employ it 
as a concrete and realistic example. 


II. MODEL SET-UP AND FORMALISM 

To set up the dark matter scenario to which our strat¬ 
egy is applied,^ we first introduce a dark sector contain¬ 
ing (at least) two DM candidates. We then assume that 
one of the DM particles is heavier than the other and the 
heavier one (henceforth denoted by Xh) communicates to 
the SM sector via the lighter one (henceforth denoted by 
X/), i.e., relevant relic abundance can be evaluated by the 
assisted freeze-out [ 18] . However, we further assume that 
X/i’s are not directly communicate to x/’s, but through 
an intermediate state A which subsequently decays into 
a xi dark sector particle a via, in general, multiple 
intermediate states. For simplicity, a single intermedi¬ 
ate state B is taken throughout this letter, but we later 
discuss the necessity of it.^ Particle a eventually decays 
into a photon pair, being taken as a source of 7 -ray ex¬ 
cess. In this sense, it could be regarded as an ALP and 
dark pion. We again refer to FIG. 1 demonstrating the 
example dark cascade scenario that we shall discuss in 
more detail. Here Xh pairs annihilate into A plus A' and 
A decays into B and a dark sector particle x'l solely for 
full generality. Since subsequent dynamics of A' and x] 


^ It is well-known that sequential decays through and into SM par¬ 
ticles typically invoke broadly-distributed energy spectra. See, 
for example, Ref. [ _] in the context of cosmic gamma-ray physics 
and Refs. [■,'1, .54] in the context of collider physics. 

^ More dedicated dark matter model building to satisfy all features 
of the generic set-up is beyond the scope of this paper. We 
instead leave it as a future work [. ] 

^ In general, xi vasty be either an unstable particle which subse¬ 
quently decays into lighter particles, or even particle a. However, 
our argument does not depend on the attributes of xi ■ In addi¬ 
tion, A and B can be generally either dark or SM sector particles, 
but we just assume that A and B are dark sector particles for 
simplicity. 
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is irrelevant to the later argument,^ we simply omit their 
further processes. 

Throughout the later argument, we take the assump¬ 
tion that particles A and B are scalars or produced in an 
wnpolarized way unless specified otherwise. With the as¬ 
sumption that Xh is non-relativistic, it is straightforward 
to have the range of the boost factor of B, xb- 


1b < lB<lt, 



where 7 a denotes the fixed boost factor of A (i.e., 
m^^/mA) and Eg {Pg) denotes the energy (momentum) 
of B measured in the rest frame oi A. As is well-known, 
the distribution in xb is flat. To develop the intuition 
on the boost of particle a, we cast its energy in terms of 
boost factor jb- 

Ea = e:xb+p: \lii-i cos e: , ( 2 ) 

where E* [P *) is the energy (momentum) of a measured 
in the rest frame of B while Ea is the corresponding en¬ 
ergy measured in the laboratory frame, and 0 * is the an¬ 
gle of its emission measured from the boost direction of 
B. Obviously, 7 q becomes 1 at cos 9* = — 1 together with 
the relation of xb = Pal'^a-, which arises in much of the 
relevant parameter space. The reason is that xb itself is 
given by a range as in Eq. (1), not a hxed value like 7 a 
so that such a special value of xb can be easily covered. 
It turns out that the boost distribution of a, g(ia) typ¬ 
ically develops a rising-and-falling shape with a plateau 
region in the middle of it as far as 7 )) is away from 1 . 
More details of g{ia) are determined by the associated 
mass spectrum of particles and underlying dynamics such 
as spin correlation if allowed. However, to avoid any un¬ 
necessary complexity later, we simply suppose that g{ia) 
takes a rising-and-falling structure starting from 7 ^ « 1 . 

When it comes to photons emitted from o, the relevant 
7 -ray energy spectrum becomes more involved. Now that 
the two photons are not in the rest frame of a, their 
energy should be Lorentz-transformed like Eq. (2): 

E^ = E;{xaPVii^cose;), (3) 

where E* {E^) is the photon energy measured in the 
rest frame of a (in the laboratory frame) and 9* is the 
angle of photon emission measured from the boost direc¬ 
tion of a. If a is a scalar, cos 9* becomes flat so that 
again the energy distribution in Ej becomes rectangular 
for a given Remarkably, E* is commonly included 
in the rectangles for any 7 ^ [52, 53]. This implies that 


one simply “stacks up” such rectangles for all possible 
values of Xa to get the energy distribution of 7 -rays, 
f{E^). This procedure is then translated into perform¬ 
ing a Lebesque-type integral over the y^’s contributing 
to the given E^ [52, 53]: 



However, in general, the above integral is not analytically 
doable, so we employ the asymptotic ansatz that was 
originally proposed in Ref. [53]. 


f{E^) ^ exp 



(E^ 


2 

Ke; 

eJ\ 


( 5 ) 


where w is the width parameter being responsible for the 
details of the distribution. Its performance and applica¬ 
bility were very successful in a wide range of examples of 
collider signatures [53, 54, 56-59],^ and therefore we an¬ 
ticipate that it will be also relevant to the case at hand. 
We are noticed that any ansatz should satisfy the sym¬ 
metry property under the operation of -O- - 5 ^, and 
thus we slightly modify the above ansatz by introduc¬ 
ing another fit parameter p determining the power of the 
argument in Eq. (5): 


fniEj) = N exp 


w 

(E, 


~2 

\e; 

eJ . 


( 6 ) 


where N is the overall normalization parameter. We shall 
conduct data ht procedures later with this modified tem¬ 
plate. 

An interesting observation should be mentioned here. 
It is straightforward to see that E* is the geometric mean 
of minimum E^ (i.e., cos 9* = —1) and maximum E^ (i.e., 
cos 0* = -I- 1) for any y^ according to Eq. (3). This implies 
that the distribution in logarithmic E^ is symmetric with 
respect to E^ = E*. In other words, the center position 
of the photon energy distribution appears as the unique 
peak in the energy distribution as mentioned before, and 
is exactly the same as E* = mal‘2., which can be used 
for mass measurement of particle a. This robust link 
between the peak position and a physical property (the 
mass of a) is rather remarkable in contrast to other DM 
interpretations of the GeV y-ray excess wherein such a 
connection is not manifest. For more generic and detailed 
discussions and physical implications about these obser¬ 
vations, see Ref. [53] and its follow-up works [54, 56-59]. 

Before moving on to our data analysis, we briefly dis¬ 
cuss “would-be” minimal scenario where particle A im¬ 
mediately decays into particle a, not through the inter- 


Particles A' and x'l could even decay into SM particles unless 
they are severely constrained by other observational data (e.g., 
excess in cosmic-ray positron measurements, [2— ■>]) or cosmolog¬ 
ical bounds. 


® Similar idea was also used for non-standard interpretations of 
the 750 GeV diphoton excess at the LHC [60]. 
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mediate state B.® Obviously, the boost distribution of 
particle a in this case is given by a box shape (see Eq. (1) 
with B replaced by a). The resulting gamma-ray spec¬ 
trum contains a plateau region in the middle of it un¬ 
less the minimum approaches 1. We find that the 
GeV excess does not show any apparent plateau struc¬ 
ture in the plane of versus dN/dEj. Of course, one 
could render close to 1, but this demands an unnatu¬ 
ral tuning among relevant mass parameters, i.e., the fixed 
'Ya{= ‘>^xh/'^A) should be close to E’^/rUa- Furthermore, 
the box-shaped probability distribution for y^ ensures a 
non-zero probability at ya = 1. This would result in a 
sharp spike at E^ = ma/‘i in the y-ray spectrum [57], 
which is not accommodated by the observational data as 
well. We also explicitly checked whether the would-be 
minimal DM scenario can explain the GO y-ray excess, 
and found that the relevant fit does not reproduce the 
data in a reasonable level. Therefore, the would-be mini¬ 
mal scenario is highly disfavored based upon this series of 
observations. Instead, the scenario delineated in FIG. 1 
can be considered minimal, while evading all the above 
unwanted situations. 


III. DATA ANALYSIS AND DISCUSSIONS 

We next conduct the fit to the spectrum of the ob¬ 
served GC GeV y-ray excess with the expected shape in 
Eq. (6). The relevant data points are taken from Ref. [17]. 
Since their data is reported in the plane of E^ versus 
E^-^, i.e., weighted energy spectra, we practically per¬ 
form the fit with E?^fM{E^). The fits are done with both 
of their regions of interest (ROIs): (i) 1° < |6| < 20° and 
|Z| < 20°, and (ii) full-sky with |6| > 1°, where b and 
I are the Galactic latitude and longitude, respectively . 
In FIG. 2, our fit results are shown for both ROIs: the 
upper panel for ROI (i) and the lower panel for ROI (ii). 
The data points are represented by black dots and their 
error bars are represented by black lines. The best-fitted 
curves are exhibited by blue curves. 

We see that our fits are in a rather good agreement 
with the observed y-ray spectrum for both ROIs. In 
other words, the data is well-reproduced by our fitting 
template. To quantify the goodness of our hts, we eval¬ 
uate the values which are 30 and 21 for ROI (i) and 
ROI (ii), respectively for 18 degrees of freedom (i.e., 22 
data points subtracted by 4 fitting parameters such as 
N, w, p, and E*) between 0.3 and 50 GeV. More quanti¬ 
tatively, our values are twice smaller than those from 
other studies such as Ref. [17]. These fit results indicate 
that our model scenario can provide a very good explana¬ 
tion for the observed y-ray excess. As advertised before, 


See Ref. [GI ] for a dedicated discussion on physical implications 
and features of this would-be minimal scenario with regard to 
the DM. 
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FIG. 2: Upper panel: the y-ray spectrum of the DM compo¬ 
nent taken from Ref. [17], corresponding to the data in ROI 
(i) (1°|5| < 20° and j/j < 20°). The ht is performed with 
22 data points, the best-fitted result is shown by the blue 
curve. Lower panel: the same thing as in the upper panel, 
corresponding to the data in ROI (ii) (full-sky with |6| > 1°). 


the fit can tell us about the mass of the particle decay¬ 
ing into two photons, m^. The best-fitted nia’s extracted 
from E* values for both ROIs are given by 

_ r 0.90 ± 0.14 GeV for ROI (i) 

“ t 0.90 ± 0.18 GeV for ROI (ii), ^ 

where the reported errors here are la statistical uncer¬ 
tainty. Note that the peak positions in the distributions 
in FIG. 2 are slightly away from mo/2 unlike the previ¬ 
ous theoretic argument because they are weighted energy 
spectra so that the peak position becomes shifted to the 
higher energy region. The fitted masses for both ROIs 
are in a great agreement to each other, more convincing 
ourselves that the data in both signal regions comes from 


a common source. 
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IV. CONCLUSIONS 

In conclusion, we have shown that for the GC GeV 
7 -ray excess, an alternative avenue of generating contin¬ 
uum bump is available with the DM interpretation. More 
specifically speaking, as a basic and minimal set-up, we 
introduced a dark sector in which there exist (at least) 
two DM candidates, the heavier and the lighter ones. 
The heavier DM particle is set to communicate to the 
SM sector via the lighter DM particle, but heavier DM 
particles are assumed to annihilate into an intermediate 
state rather than directly into lighter DM particles. This 
intermediate state further decays into a lighter DM par¬ 
ticle together with an ALP or dark pion, which directly 
decays into a pair of photons, via intermediate dark sec¬ 
tor state(s). In this scenario, we showed that the first 
intermediate state comes with a fixed boost factor, and 
in turn, the second intermediate state comes with a flat 
boost distribution. Thus, the ALP or dark pion gets 
varying boost factors characterized by a rising-and-falling 
structure, naturally resulting in a continuum photon en¬ 
ergy spectrum even though it decays directly into two 
photons. 

One of the novel features in this DM interpretation 
is that the relevant energy spectrum has a symmetry 
property with respect to the rest-frame energy of pho¬ 
tons (A*), which motivates and constrains the form of 
the ansatz describing the relevant energy distribution. 
Therefore, the mass of the particle decaying into a pho¬ 
ton pair can be measured by reading off E*, which is 
given by half the mass of the ALP or dark pion, from the 
fit procedure. In addition, considering the DM scenario 
to be employed here, one could take the GeV excess as 


a “smoking-gun” signal for the non-minimal dark sector, 
which currently receives increasing attentions. 

With the given model set-up, we have introduced a 
template for the GG GeV 7 -ray excess. We then have 
studied how well such a template can reproduce the ob¬ 
served energy spectrum of GG 7 -rays. Performing the 
standard fitting procedure, we have obtained = 30 
and 21 for ROI (i) and (ii), respectively, and both ROIs 
have consistently reported that the mass of the ALP or 
dark pion is ~ 900 MeV. We expect that these values will 
be better-established as the Fermi-LAT collaboration ac¬ 
cumulates more cosmic-ray data. 

In forthcoming work [55], we shall construct realistic 
DM models to accommodate the basic set-up that we dis¬ 
cussed in this letter. In particular, we shall investigate 
phenomenological and cosmological implications of the 
models in the aspects of direct detection, collider phe¬ 
nomena, and dark matter relic density. In addition, we 
shall study the connection between underlying DM model 
parameters and the fit parameters other than E* in the 
context of a concrete DM model. 
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